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SECTION I 


Introduction 


Previous theoretical [1-3] and experimental studies [4] established 
the foundation for a correlation between the plane strain fracture 
toughness k lc and the ultrasonic factors v L e s /m through the 
interaction of waves with material microstructures, i.e. grain size or 
second-phase particle spacing depending upon the material system. This 
suggests a means for ultrasonic evaluation of plane strain fracture 
toughness. Ultrasonic methods can therefore be used not only in the 
evaluation of material properties such as moduli and porosity but also 
in the fracture properties such as crack size and fracture toughness. 

There are situations where material properties keep changing. 

Stress induced deformation during manufacturing processes often cause 
continued changes in material properties. It is well-known that pre- 
cipitation, martensitic transformation, void nucleation, etc. often lead 
to changes in material properties. Fracture toughness enhancement has 
been observed in a number of ceramic systems due to stress induced 
transformations e.g. Zirconia (Z r 0 z ) particles contained in ceramic 
matrix are observed to transform from a tetragonal to a monoclinic 
crystal structure at sufficiently high stress environment. Circumferen- 
tial microcracking may occur and a process zone is usually formed in the 
vicinity of the stress raiser to reduce crack-tip stress intensity. 

Precipitates, martensites and dislocations are typical examples of 
inhomogeneous inclusions [5,6,27,32], i.e., regions with distributed 
transformation strains alias eigenstrains. Since the scattered 
displacements are directly related to the eigenstrains [7,8], the 
scattered field measured can be used to studied the changes caused by 
the presence of precipitates, martensites, voids, etc. 
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The purpose of this work is to lay a foundation for the ultrasonic 
evaluation of microcracking. The work establishes an "average" theorem 
and wave scattering effects of a microcrack and disbonded microvoid. 
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SECTION II 


Technical Background 

Definition of Microcrack Nucleation 

A microcrack represents solid-vapor surfaces in a material and it 
has a dimensional scale comparable to some microstructural features, 
e.g. the grain size. The use of the term microcrack should not be con- 
strued to imply brittle fracture. An microcracking event is obviously 
the origin of fracture. In what follows, a brief review of different 
mechanisms of microcracking is outlined. 

Microcracking in Polycrystals 

In single phase crystals, a microcrack nucleated from defect loca- 
tions can grow at the competition between dislocation generation and 
bond rupture of the crack tip. The activation or generation of disloca- 
tions that pile up at the crack tip essentially determines whether the 
crystal is brittle or ductile. 

For cleavage-prone polycrystalline materials with brittle parti- 
cles, such as A533B steel, the brittle fracture proceeds by a mechanism 
involving a slip band blocked by a carbide. The most heavily strained 
carbide particles are usually found to be a few grains ahead of the 
crack-tip. The microcracking event therefore occurs at that location 
under suitable loading conditions. For a polycrystall ine material with 
a matrix that resists cleavage the microcracking can occur by void 
nucleation and growth mechanism. Stresses arising due to the incompati- 
bility between the inhomogeneities and the material matrix may cause 
eventual cracking of the particle (inhomogeneity) or the interface. 

Depending upon the "size" of the particles involved, secondary 
plastic zone near the inhomogeneities may or may not be of importance. 

In essence, the continuum plasticity models can be used for large 
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i nhomogeneities, 6^1 ym detail dislocation structures must be 
considered for the medium size particles, 0.01 ym < 6 < 1 ym and at the 
decohesion of interface, the size of the plastic zone seems to be de- 
pendent on the diameter of particles rather than the strains. Shearing 
of the weak inhomogeneities may be an important mode of microcracking in 
the case of very small particles, <5 > 0.01 ym. It is clear from the 
discussion above that the microcrack nucleation mechanism involves the 
situations of either (1) the creation of new surfaces by tension or by 
shear at a weak second-phase particle or (2) the nucleation and growth 
of voids at a strong second-phase particle. 

Microcracking in Ceramics and High Temperature Materials 

The microcrack nucleation is an event of crack initiation. The 
presence of a single microcrack is usually not of serious consequences. 
It may, however, grow in size to become a macrocrack. Multiple 
microcracks may nucleate near the vicinity of a macrocrack to form a 
microfracture process zone thus enhance the coalescence of microcrack 
that add to the length of a macrocrack. Voids nucleated along grain 
boundaries may be flatened out into cracks under favorable conditions. 
Spacing between the voids, can be strongly dependent both on distance 
from the crack and on the duration of loading [32]. One form of 
transformation toughening of ceramics with intercrystall ine ZrO£ dis- 
persions was found to be largely caused by microcrack nucleation and 
extension. The phenomenon of fracture toughness enhancement was ex- 
plained by the use of circumferential microcracking model [30] and by a 
continuum elastic plastic fracture mechanics approach. 

Dynamic Models of Microcracking 

A simple dynamic model for microcracking was first proposed by Vary 
[4] where he considered the stress wave energy required to create a 
microcrack of the diameter same as the second-phase particle. He 
considered this potrion of the energy to be totally consumed by the 
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particle during its fracture and gave an interelation between the 
ultrasonic and fracture toughness factors. A theoretical verification 
of the model reequires information about the scattering of an embedded 
flat microcrack. In the fracture of ceramic materials, there are cases 
where the failure prediction requires the direct evaluation of the 
maximum stress intensity factor by either proof testing and/or acoustic 
emission. 

For the elastic wave scattering from planar cracks was given a 
starting point for the integral formulation is the representation of the 
scattered fields in terms of the crack-opening displacement, auj, 

u i> S> = - C jkrs f Au d n k dS' (1) 

*S 

Numerous investigators have examined the quasi-static scattering from 
elliptical and circular planar cracks. Since the crack-opening 
displacement, AUj, is not a known-priori in Eq. (1), approximations 
can be obtained by using well-known exact solutions for cracks under 
static loading. The obtained solution are thus appropriate for long 
wavelength. An estimation procedure for fracture crack parameters was 
previously given by Budiansky and Rice [11]. On the other hand, Fu, Co 
and Dzeng [14] presented a theory for elastic wave scattering and crack 
sizing which is appropriate at wider frequency range. Both of these 
studies were in the frequency domain. It should be noted that informa- 
tion in the frequency domain and the signals in time domain are related 
via Fourier transforms. The solution presented here applies to the 
frequency region corresponding to dimensionless wavenumber ka from zero 
to about less than 2 . 

Plan of Study 

The microfracture problems can be analyzed by first computing the 
appropriate microstress distributions, then estimating the stress 
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Intensity factor K for typical mlcrostructural defects located In the 
stress field, and finally equating K to the local material resistance to 
microfracture. The microfracture problem of greatest practical Interest 
concerns the grain shape, size, and relative orientation that dictates the 
onset condition for boundary microfracture between grains. 

As a beginning In this program, the elastic wave scattering from a 
planar crack under Incident plane waves Is studied. The specific 
geometries studied are the penny-shaped circular crack and the elliptical 
crack. Far-fleld quantities such as the reduced scattering amplitude that 
are related to the shape, size and orientation of the crack will be 
determined. Near field quantities such as the crack size, the stress 
Intensity factors and the crack opening displacements are briefly studied 
[13,14]. The approach taken Is based upon the methods developed in 
[7,8,12]. 

As a second part of the work for this program, an average theorem Is 
developed to study the dynamic effective properties of media with randomly 
distributed Inhomogeneities. The approach taken Is deterministic In 
nature and the solution Is based upon the elastic wave scattering of a 
single Inhomogeneity, a self-consistency scheme and an averaging theorem 
over strain energy and kinetic energy [8,15]. 
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SECTION III 


Theoretical Studies 

The aim in non-destructive evaluation (NDE) by sound and ultrasound 
is to extract material and fracture properties from pulses sent and re- 
ceived at transducers. The analysis is usually given either in the time 
domain or the frequency domain. The introduction of the size factor, 
the wavelength, in the model mechanics problems in NDE, plays the impor- 
tant role of relating the far-field measurable quantities such as phase 
velocities and attenuation to the near-field physical situation such as 
local geometric dimensions and elastic properties. This aspect of the 
analysis thus allows "nondestructive" testing methods to be employed for 
"experimental" validation of the predictions by the theories. 

Solution to the elastic wave scattering due to a single embedded 
inhomogeneity is available by different methods that are appropriate at 
different frequency ranges [16,10]. The methods that offer a solution 
in an analytic form and are useful for inhomogeneous media with multiple 
components are the longwave approximation, the polarization approach and 
the extended method of equivalent inclusion [7,8]. 

There are several averaging schemes or theorems that exist in the 
literature for finding the dynamic effective moduli and mass density ._ 
While some efforts concentrate on the average stress a and strain e 
fields or on the average displacement field u, others concentrate on a 
variational approach. These theories are appropriate mostly at Rayleigh 
or longwave limits and do not exhibit dispersive effects. Dispersive- 
ness and attenuation are important in dynamic material properties 
evaluation. 

In what follows, the scatter of elastic waves due to a thin flat 
ellipsoidal inhomogeneity is first studied. Penny-shaped and ellipti- 
cal ly-shaped cracks are studied as special cases. Fracture crack 
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parameters are presented. An average theorem appropriate for dynamic 
effective mass density and effective moduli is developed via a self- 
consistent scheme. Effective material properties of two-component media 
consisted of randomly distributed spheres are given as a special case. 


Prel iminaries 

(1) Displacement field due to the presence of mismatch 


The inhomogeneous media considered in this paper are assumed to be 

consisted of a homogeneous matrix, of elastic moduli C° and mass 

density p°, and a distribution of inhomogeneities with moduli and 

(r ) 

mass density p occupying regions ^ r , r = 1,2, n, Fig. 4. The 

total displacement field u can be separated into two parts, 


u 




( 2 ) 


where the superscripts (i) and (m) denote "incident" and "mismatch," 
respectively. It is clear that when no "mismatch" components, i.e. 
inhomogeneities, are present, the total displacement field is entirely 
the same as the "incident" displacement wave field. On the other hand, 
if there is no incident wave field, the only field that exists is the 
null field. 

(2) Eigenstrains 

The eigenstrains e* which are also termed transformation strains 
are defined as the part of the total strain e that must be subtracted 
before the remaining part can be related to stresses g through Hooke's 
law, 
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e - r" 1 
e rs L jkrs a jk 


( 4 ) 


where C - * are the elastic compliances. The method of equivalent 
inclusion is a method that allows the inhomogeneous media be replaced by 
media of homogeneous matrix effective moduli with distributed trans- 
formation strains in the regions originally occupied by inhomogeneities, 
hence 


1 0 i n n 

e*s in fi r , r = 1,2. ..n. 


(5) 


For the two prolems to be equivalent, the transformation strains must 
give rise to a field that is exactly the same as that of the "mismatch" 
field, u^ m K This leads to equivalence conditions that insure 
identical field quantities at any given point be obtained in the two 
problems" 1 ' [9]: 


A ^jkrs 


u “< f)+C Jkrs ^s 1)(F) = - 4C < 


jkrs 




in Q 


( 6 ) 


Ap to 


u! m) (F) 


C jkrs 


= 


'rs 


^(r) = -Ap 


(0 


Uj 1 ^ (r) 


in Q 


(7) 


There are two types of transformation strains or eigenstrains that 
arise in el astodynamic situations due to the mismatch in elastic moduli 
ac and mass density Ap . It is often convenient and useful to define 
associated quantities such as 


m jk C j krs 



( 8 ) 


"•"The conditions (6.7) are similar to those of Willis (1980) and 
those of Mura, Proc. Int. Conf. on Mechanical Behavior of Materials , 5_, 
Society of Materials Science, Japan, 12-18 (19/7). 
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*. = c e *(2) 

J L jkrs rs,k 


( 9 ) 


where mj|< and tt j are referred to as moment density tensor and 
equivalent force or eigenforce, respectively. 

(3) Volume average and time average 

The volume average and time average of a field quantity say F(r,t) 
is denoted by using brackets < >, and < >j, respectively, and are 
defined as 


<F(r,t)> = 1 / F(v,t)dV 


( 10 ) 


<F(r,t)> T = | / F(r,t)dT (u) 

where V and T stand for volume and time period, respectively. 

(4) Volume integrals of an ellipsoid associated with inhomogeneous 
Helmholtz equation 

Volume integrals of an ellipsoid associated with the integration of 
the inhomogeneous Helmholtz equation are used in this work. The 
inhomogeneous scalar Helmholtz equation takes the form: 


2 2 

v o + k $ = -4 tty(j2) 


( 12 ) 


where Y(jr) is the source distribution or density function, v 2 and k are 
the Laplacian and wavenumber, respectively. A particular solution to 
Eq . ( 12 ) is 


$(r) = / y(jr' )R"^exp(ikR)dV' , R = | 

a 


(13) 
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in which (4irR)-lexp( ikR) is the steady state scalar wave Green's 
function and n is the region where the source is distributed. The 
source distribution function y{r) can be expanded in basis functions or 
polynomial form, depending upon the geometry of volumetric region fi. 

For an ellipsoidal region, the choice of using a polynomial expansion 
separates this work from other theories of elastic wave scattering: 

Y(r') = (x')V) 1 'U , ) v (14) 


in which x, y, v are integers. 

For elastic wave scattering in an isotropic elastic matrix, two 
types of volume integrals and their derivatives need to be evaluated: 

ip (jr) = f R -1 exp(iaR)dV‘ (15a) 

n 

* k (r) = / x k R _1 exp(iaR)dV' (15b) 


*k£...s^ = ' x k x r 


,X s R 


exp(ioR)dV* 


(15c) 




(16a) 


♦k.pto (16b) 

• • • 

where a 2 = w 2 /V 2 = pto 2 /(x + 2p). The other type, the ^-integral s, are 
obtained by replacing a with B in Eqs. (13,14), where B = pw/(my) 
Details of the integration are given in [12]. In this work only 
limiting value at r-K) and r-*» are of interest. 
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An Isolated Flat Ellipsoidal Crack 
(1) Formulation [28,13,8] 


Consider the physical problem of an isolated inhomogeneity embedded 
in an infinite elastic solid which is subjected to a plane time-harmonic 
incident wave field as depicted in Fig. 1. Replacing the inhomogeneity 
with the same material as that of the surrounding medium, with moduli 
c jkrs and mass density p, and include in this region a distribution of 
eigenstrains and eigenforces, the physical problem is now replaced by 
the equivalent inclusion problem. 

The total field is now obtained as the superposition of the inci- 
dent field and the field induced by the presence of the mis-matches in 
moduli and in mass density written in terms of eigenstrains and 

eigenforces, 7 t* ^ 

J 


F = F (i ) + F (m) 


(17) 


where F denotes either the displacement field uj, the strain field 
e ij, or the stress field o^. The superscripts (i) and (m) denote 
"incident" and "mis-match", respectively. 

For uniform distributions of eigenstrains and eigenforces, the 
fields can be obtained as: 




U m ( r ') = - 

m — 


$•>') - C, 


, 0 ) 


J jnr_ ■ 


jkrs rs 


S jm,k'(l) 


(18) 




+ u 


(m) 

n,m 


)/2 


(19) 


o (m) = C 

pq pqmn mn (20) 
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where a comma denotes partial differentiation and 




’a 9jn, ( ^' )dV ' 


( 21 ) 


in which gj m is the spatial part of the steady state elastic wave 
Green's function and q is the region occupied by the inhomogeneity. It 
is noted that the integrals Sj m and their derivatives must be 
evaluated for the region _r > _r‘ and for the region r < r' , Ref. 12. 

The solution form represented in Eqs. (15-19) gives the fields inside 
and outside an isolated inhomogeneity of arbitrary shape. 


(2) Far-field scattered quantities 

Let the incident displacement field be longitudinal and of 
frequency w , amplitude u 0 : 

U J'^ = U o q j ex P[i aX -jk-j - Iwt] 


( 22 ) 


where i^ = -1 and qj is the unit vector in the normal direction of 
the plane time harmonic wave and k-j is wave vector. For a linear iso- 
tropic medium, the spatial part of the free space Green's function is 
well known. Substituting gj m (r-r‘ ) in Eqs. (19,21) and using the 
limiting concept 13: 


limit a_ if* - A. constants 

a 3 +0 J J J 


(23) 


limit a_ = B. . constants 

a 3 ->0 J 1J 1J 


(24) 
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the scattered displacement i4 s ^(r,t) from a thin elliptical flat crack 
can easily be obtained as: 


( “ a l> 3u o 


(c, a] )V 


r oo 


(25) 


= [(C G m expi ar)/ar + (D H m expi er)/er]exp(-icot) 
where 

G m = - (a 2 /a l>[-V/j + (1 - 2c ‘ 2/B2 ) t m B jj +2 (“ 2 / |s2 ) t m ,! k i j B kj :i 1261 

H m " (^A,)[-(B/«) 3 (Vj- { m j) A j 
- 2(8/»)\B* m + 2(B/a) 2 Vk£j B^.] 
c = [j 0 (ar e )+j 2 (ar e )]/3 = sinar e /(ar e ) 3 - cosar e /(ar e ) 2 (27) 

D = [J 0 (6r e )+j 2 (er e )]/3 = si nBr e /(er e ) 3 - coser e /(er e ) 2 


in which m,j,k = 1,2,3, and (ai,a 2 ), “>2, denote the 

semi-axes of the flat ellipsoidal crack, direction cosines of scattered 

displacements, longitudinal wavenumber, shear wavenumber, respectively. 

Also, Aj and B j ^ are the reduced non-dimensional form of Aj and 

Bj|<, respectively, defined as follows: 

A j = Aj/(Apu 2 u 0 ) , Ap = p' -p , (p‘ = 0 for a crack) (28) 
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(29) 


Jk = -B jk /(iau 0 ) 

Expressions for the differential cross section dPM/dn and total 


cross section P(w) can be obtained as [19]: 

a L (e ,4>) + (a/3)c T (8,i,) (30) 

P(u) = / z [o L (e,<j)) + (a/f5)a T (e,<j>)]dft (31) 

where dn is the differential element of solid angle and 

a^a^(e,<J>) = (cta 1 ) 6 [C G m ][C G m J (32) 

B 2 a T (e,<fi) = (aa 1 ) 6 [D H m ][D H m ] (33) 


in which the super bars denote complex conjugate. It is noted that the 
constants Aj and Bj k must be evaluated from the equivalence 
conditions Eqs. (6,7) with the use of the limiting concepts in Eqs. 
(22,23) and of the integration method developed in Ref. C12D. 

(3) Determination of Aj and 

In Eqs. (22-26) the scattered displacement field is given in terms 
of the "reduced" form of the eigneforces and eigenstrains , i.e. Ajf and 
Bjfc. These constants must in turn be determined from the equivalence 
conditions. Writing the incident wave field in a Taylor series the gov- 
erning simultaneous algebraic equations can be easily obtained. Since 
f -j j [0] and F mi -j[0] vanish automatically, these governing equations 
become uncoupled and lead to a three by three system for Aj and a six 
by six system for Bj k . For a linear elastic medium, they are: 
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(34) 


*A.V ma 1 + A s = - q s 

{4X5 stW° ] * 24pD stjk C0 ™jk + < X5 S t *L + 

= - tAX{ 5tVm + 24uC l s V 


where the subscripts s,t,m,j,k = 1,2,3 repeated subscripts denote sum 
from 1 to 3, and 


4ir PM 2 f js (r) = - 3 2 4> 6j s + 


mj 


) * • 
mj 


(36) 


4 " p “ 2 D stjkW - 2 " c *-syit + ♦•stjk 3 

- + p . js « kt 3 (37) 

' x “ 2 *'mn 6 jk 

in which 

AX = -X, Ap = -p, Ap = -p 


The <|>- and ^-integrals and their derivatives are evaluated by the 
method suggested in Ref. [121 Retaining terms up to (aai) or (fa^) 
of the fourth order, the constants are obtained as: 


A* = - q./{u o « 2 ApfjC03 + 1} 


no sum on j 


(38) 


f.ioi = (f^ cod, f 22 co:, f 33 co:) 


(39) 
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(40) 


{B}} = Cb ij ]" 1 Cc j : 1,j = 1,2,3 

{B*.} = - q .q ./Cl + t(<P- CO] + <i.COD)D, no sum on i,j. 

1 J 1 J J 

i t J 

1 , j = 1 ,2,3 5 = 1/4 tt 

where in Eq. (37) = B^ » B^ = B|^ } 

Cj = (AX + 2Ay)C1 + q^D > J — 1 >2,3, 

b-|-| = (x + 2p) + Xc^ » j j C03 + 2(x + 2p)?({) , -j -j COD + 2p^<j/,-|iC0D 

b^ = X + Xc^jjjCOD + 2pCif) ,|-| C0D+ 2£Xc{>, 22^03 

b-^ = X + X?^,jjC0D + 2p5'1' >11 COD + 2XC4> , 33 ^-®-* 

b£i = X + X?i|;,jjC0D + 2y5^, 22^03 + 25X<j> ,-| -j COD 

b 22 = U + 2u) + Xtf>,jjC0D + 2p^, 22 C0D + 2?(X + 2y)<f> , 2 2 Co:! 

bg^ = X + X^ip , j jCOD + 2p?^, 22^03 + 2?X4> j^gCOD 

^31 = ^ + XC^>jjCOD + 2y?^,23C0D + 2X?<i> , 1 -| COD 


(41) 

(42a) 

(42b) 

(42c) 

(42d) 

(42e) 

(42f) 

(42g) 

(42h) 
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b 32 = X + X^-jCOH + 2 u ^, 33 C0] + 2\^, ZZ 101 


(421 ) 


^33 (x 2y ) + XW.jjCO: 2y^^» 33 L0] + 2 (x 2ij)54> » 3 3 C03 (42j) 

Note that . f b^. . In Eq. (38), ^ <J> 2 = <P Z 2 C0], cj> 3 = <J> 33 E0D. 


The f- , 4>- functions are given as: 

4-rrpu)^ fj s [0] = “ Picons j s + t|»,j s t03 “ ^’js 1103 (43) 

AcOH = ua 1 a 2 e Z {I 0 - ((Ba 1 ) 2 /16)I 1 + i(8/3)p} (44) 

cf> all C03 = -C(Tra^a 2 B 4 )/12)3. (45) 

4», 22 C0]= -C( 1 ra 1 a^B 4 )/12]. (46) 

<t>, 33 E0: = 0 (47) 

i n which 



diji 

fiW 


F(0,k). f- 
a l 


(48) 


y°° ^difj _ 2_ j- E(e,k) 

0 (a 2 + *)a(*) ' 3 1 k 2 


k |2 

- V F(e,k)} 
k 


(49) 
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(50) 


I 



_ 2_ |- F(9 ,k) 

(ag + ^)A(ip) a l k 2 


E(e,k) 1 

,2 } 


F = / 0 o d(iJ o 1/? , E = / 0 (l-k 2 sin 2 a)) 1/2 da) (51) 

0 (l-k 2 sin w) ' 2 0 

and as , 9 ^tt/ 2 , k -*■( 1 -a 2 / a-j ) and k — ( 1 -k )- >- a 2 / a-j , if ^ ^ 2 . 

a^O and a-j = a 2 , we have I q = ir/a^ ^ = I 2 = u/Za^ The 4^ - functions 
are obtained by replacing 9 with a in the $ - functions. 

(4) Stress Intensity Factor (SIF) and Crack Opening Displacement (COD) - 
preliminary results. 

The general expression for the total displacement field can be 
easily obtained, by using Eq. (2) and the expression for the mis-match 
displacement field [9], as follows: 


u m (n) ■ \,,Cr)* j ♦ W r,A jk+ 

+ ^mij^— + ^mijk^— ^ijk+ — 


(52) 


To find the crack face opening displacement the functions f K and F< 
must be evaluated for r_ inside the crack. 

In this section, the solution for a crack is that degenerated from 
an ellipsoidal void by letting a3+0. The crack opening displacement 
for an elliptical crack can be defined as 


Au. = 2u . = h(ei. + 77 *) (53) 

1 1 ol 1 

where h is the half thickness of the ellipsoidal void in the 
X 3 ~di recti on: 
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(54) 


h = a 3 (l - x 2 /a 2 - x 2 /a 2 ) 1/2 


For the case of uniformly diistributed eigenstrains and eigenforces, 
the crack opening displacement is, therefore, obtained as: 


or 



+ B 3j )a 3 ha a 


3+0 




1 

.2 


4 $ 
4 


A 


Ap 

- A ? + B 5i 

P 


(55) 


where u 0 is the amplitude of the incident wave field and Aj and 
B 3 j are given by Eqs. (38-41). 

The stress intensity factors can be easily obtained by finding the 
strain energy change W, i.e. the interaction energy, due to the 
presence of the crack and then differentiate with respect to the crack 
dimension. An approximation for low frequency range can be obtained by 
using the SIF-COD relation [11]. Surface-wave interaction with an edge 
crack was studied in [33]. 

(5) Numerical calculations and graphical displays 

It is clear from Eqs. (38-41) that the uniformly distributed eigen- 
strains and eigenforces, in their reduced form, Bjj and A j , respect- 
ively, are dependent only upon the incident wave field characteristics 
and the geometric factors of the inhomogeneity for a given material 
system. It is observed from Eqs. (38-42) that Bj and Bjj as 
functions of dimensionless wavenumber of a^ or aj would exhibit 
large peak values at certain incident wave frequencies when bjj 0 . 

The values of these critical frequencies depend only upon the matrix 
elastic moduli, the crack dimensions and the measurement direction. For 
a given crack aspect ratio, ^ 2 /^ 1 , and a measurement direction, 
the difference in frequencies at subsequent peak values is proportional 
to a^, the largest dimension of the crack. Fig. 9. 
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Computational data of elastic wave scattering due to a flat crack 
embedded in any given isotropic material system can be obtained by 
employing Eqs. (25-29) and (34-51) with Ap = -p, AX = -X and Ay = -y . 
Scattered displacement amplitudes can easily be obtained for any given 
crack aspect ratio. 

(a) Elliptical Crack: In Figs. 7-9 computational data for back 

scattered situation are displayed for a crack in aluminum of different 
crack aspect ratios, i.e. a 2 /ai = 1, 1/2, 1/10 at various wave 
incident angle. Critical frequencies are observed to be present in 
B*j but no critical frequencies are identified for Aj, j = 1,2,3 
Figs. 10-13. The critical frequencies for C|G m | and D | H m [ are, 
however, clearly identified. The position along the aa^ axis at which 
the first critical frequency occurs depend upon the crack aspect ratio 
and matrix elastic moduli. Sufficiently small increment in aa^ must 

be used in order not to miss any peak values. Since the soulution form 
given in Eqs. (25-27,38-51) is analytic in frequency, this can easily be 
achieved. It is noted that most of the scattered energy is carried by 
the transverse components of the scattered displacement, i.e. 

|DH m | > |CG m |, Figs. 14,15. 

(b) Penny-shaped Crack: In Fig. (7), the reduced longitudinal back 

scattered displacements at 0 = <f> = 45° and the same crack aspect 
ratio, a 2 /ai = 1, are shown for a crack in structural steel. The 
change in the absolute value of the crack size causes a compaction of 
the peaks in j C G m J along the axis of dimensionless wavenumber aaj_. 

It is of interest to observe that the difference in frequencies at 
subsequent peak values, e.g. the first and the second peaks, Afj_, Is 
inversely proportional to a certain power of a^. This fact can easily 
be algebraically demonstrated by using Eq. (42) to find the critical 
frequencies, i.e. finding the roots of the following algebraic equation: 
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(56) 



+ C£W 


= c 


3 


where the constants Ci, C 2 and C 3 are obtained by directly 
substituting Eqs. (42-51) into Eqs. (38,41) for each component of bjj. 
Since 

Af = A(aa 1 )v L /a 1 ( 57 ) 

a scale in frequency is also given in Figs. (10-20). It is clearly seen 
from these figures that for small cracks, accuracy of solution at lower 
dimensionless wavenumber is advisable and for larger cracks a wider 
validity along the axis of dimensionless wavenumber is of interest. It 
should be noted that the peaks in these figures are very large. 

In Figs. 14, 19 and 20, computational data for back-scattered 
reduced transverse displacement amplitude are displayed for different 
aspect ratios, a 2 /ai = 0.01 and a 2 /ai = 0.50. For a fixed 
aspect ratio, again, it is clear that the difference in frequencies at 
subsequent peaks, say Afj or f^, is inversely proportional to the 
third power of a^, the largest dimension of the crack. All calcula- 
tions are done for given aaj. Note that the absolute scales in 
frequency are quite different for different a^ with same aspect ratio. 

Dynamic Moduli and Damage of Composites 

Consider the problem of the inhomogeneous media, as illustrated in 
Fig. 4, under a plane time-harmonic indent wave field, Eq. (22). The 
true composite thus occupies the whole region and possesses effective 
moduli C* and mass density p*. To determine the effective moduli and 
mass density, an average theorem over strain energy and kinetic energy 
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is used. The effective properties are found to depend upon a fourth 
rank tensor A and a second rank tensor D. A self-consistent scheme is 
then developed for the determination of these tensors. 

(1) Average theorem 

To determine the effective moduli and mass density, the following 
definitions are used: 


<a> = C*<e> 

<cre> = C*<e><£> 
<pV> = p*<v> 

cpv_*v^> = p *<v><v> 


(58) 

(59) 

(60) 
(61) 


where g, e, _v are the stress, strain and velocity fields, the underlined 
- denotes a tensorial quantity, and the angled brackets <> denote the 
volume average of a field quantity, Eq. (10). The left-hand-side and 
the right-hand-side in Eqs. (59) and (61) can be shown to be equivalent 
under the so-called Hill's condition E201, <oxe> = <oe>. From Eqs. 

(60) and (61) it is clear that the kinetic energy per unit volume of the 
effective medium can be made equal to that of the physical medium if and 
only if a frequency dependent mass density is defined. It can also be 
seen, that this is the same as requiring by again using the Hill's 
conditions <p,_v*v> = <pv><v> the average linear momentum per unit volume 
be the same as the effective linear momentum per unit volume, Eq. (60). 
These condttvons, Eqs. (60) and (61) are not met if the static defini- 
tion of effective mass density is used. 

Let f r denote the volume fraction of the (r)th inclusion 
material, then the volume average of the stress and velocity fields <* 
and u, respectively, are 
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<CT> 


( 62 ) 


n 



(r) 


<lj> = E f 
r=o r 

where = C^ r M r ^, (r) = 0,1,2, n, 

and for time-harmonic situation 

u = v = -iu)U 


(63) 


(64) 


Equations (59) and (61) can be rearranged, by using Eqs. (62) and 
(63), as 


<ae> = C°<e><e> + E f c (r) e (r) 
r=l r ~ ' 


<£> 


(65) 


<p V* V> = p <v><v> 


+ E Op^M r ^<V> 
r=l r 


( 66 ) 


Consider now the case when the solution form possesses a linear relation 
between the velocity and strain fields in the (r)th component and the 
average velocity and strain fields of the effective medium, i.e. 


§ 


(r) 


= A< r > 


<e> 


(67) 



( 68 ) 
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where and are tensors of fourth and second rank that must 
be determined with a suitable scheme. The substitution of Eqs. (67) and 
(68) in Eqs. (65) and (66) leads to 


<ge> = {C° 


+ E f r AC (r) A (r) }< e x e > 
r=l r 


(69) 


<pvv> = {p° 


+ 


E f Ap^V r ^}<V><V> 
r=l r ~ 


(70) 


A comparison of the above with Eqs. (59) and (61), the expressions for 
the dynamic effective moduli and mass density are obtained as follows: 


C* = C° + E f 
~ “ r=l r ' 


(71) 


P* = p 0 + E f Ap (r V r ) = dW/3 (72) 

p = 1 * JJ 

It is important to note that the tensors A^ r ^ and p( r ) are frequency 
dependent and they replace the static expressions when frequency 
approaches zero and when proper cares are taken. The assumption of a 
general linear dependence between vj r ) and <v>, Eq. (68), will have to 
be specialized such that the second rank tensor D will degenerate into a 
scalar. This specialization is automatic for randomly distributed 
spheres where = D<5 m j, in which D = Du = D 2 2 = 033* 

To determine the explict form of the tensors A^ and p( r ), the 
strain and velocity fields in the (r)th component are determined by 
using the method of equivalent inclusion as presented in the previous 
section. 
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(2) A self consistent scheme for the determination of the effective 
properties 

The tensor fields A (r ^ and are of ranks four and two, re- 
spectively, and they are functions of wavenumbers, geometric properties, 
and effective and inclusion material properties. Let the average strain 
and velocity be the same as those derived from the incident wave field, 
then the governing conditions for determining A^ r ^ and D( r ) for the 
(r)th inclusion are simply obtained by rewriting the equivalence condi- 
tions for an effective medium with mass density p* and moduli C*: 

c (r) (e (, V" ) ) .= c*(s (1) + 6 (m) -s*) , in (73) 


p (r l 2 (u 


(i) +y (m)) = 


p*u^(u ^ 1 ) - j, in Q 


(74) 


where the superscripts (i) and (m) stand for "incident" and "mis-match", 

respectively. Slightly different approaches for finding effective 

moduli for hetrogeneous material are given in [20,23] that apply to 

static cases with different constituents and situations. The strength 

and fracture responses of composites have been investigated Q7 18 
20-233 . ' - ’*■’ 


(3) Effective properties of two-component media: randomly distributed 

spheres 

Let the moduli, mass density for the matrix, inclusion material 
particles and the effective medium by denoted by Cj krs , P o ; Cj krs , 
p, » an d Cjkrs and p*» respectively. Using the elastodynamic solution 
for a single ellipsoidal inhomogeneity, the displacement and strain 
fields inside an inhomogeneity are found, when the average displacements 
are made equal to the incident time-harmonic plane wave field, to be: 
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(75) 


Eu (m) (r)] = D(r)u (i) - 5(r)u (i) , in 

Ce^ m ^(£)] = A(r)e^ - 6 ( r_) u ^ ^ , in (76) 

in which the superscripts (m) and (i) denote "mis-match" and "incident", 
respectively. Employing the volume averaging process as described by 
Eq. (10) and substituting in Eqs. (65,66), the effective properties are 
easily defined as follows: 

P * = p + fApD, D = h.. (77) 

^ J J 


C* = C + fACA 


(78) 


where f is volume fraction of inclusion material. The tensor fields D 
and A are: 

°mj = - <f mj^ )>/{f MJ [:0:i + 4tt^ p ' -p*)iu 2 } , no sum on M,J. (79) 


i n 


\,npq ■ <F radk,n^> + F ndk,mt )>S J>W 2p *“' 
which the tensors f and F are defined as: 


(80) 


4*Pw 2 f mj (r) =-As mj 


+ 4 , 


5 mj 


- (pi 


mj 


(81) 


4iTpa)' 




6 . . 

m ij 


+ 2yS* 


1 ’ i 5 mj 


2vip 


’mj 


+ ^>’mj J 


The <f> - and -integral s given in Eqs. (15,16), etc. are the volume 
integrals associated with the inhomogeneous Helmholtz equation. They 
can be carried out for an ellipsoidal region by expanding (exp ikR)/R in 
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Taylor series expansions with respect to r' , for r > r' and with respect 

to r for r < r'. Here k can be either a or 6. Details are given in 

Ref. [121. This type of expansion for the integrand is particularly 

useful in determining the coefficients of a "polynomial" distribution of 
* *( 1 ) 

77 j and £-jj . 

The fourth rank tensor Sj^pq is the connecting tensor between the 
eigenstrains and the applied strains, i.e. 


e *0) _ c (a) 
jk jkpq pq ’ 


S jkpq S jkpqV a i ,( ^ 


(83) 


for the case of uniform eigenstrains and eigenforces. In developing 
these expressions the volume average of the <Mntegral s , must be 
evaluated. Finally, it should be noted that p* and C* are complex where 
the real and imaginary parts are associated with the velocity and 
attentuation, respectively. 


(4) Example: spherical inclusion materials 

Let the spherical inclusion materials of radius "a" be randomly 
distributed over the whole volume of the matrix. If the matrix and the 
inhomogeneities are isotropy, the effective medium is also isotropic. 

It is straight forward to show that 

D mj = Vi { " < f 33 ( I )>/f 33 C0:i + 47r(p ’ ' p * )a,2] (84 

V D 

and 


S jkpq S kjpq s jkqp ~ S pqjk 
S lll S 2222 = S 3333 = C 1 
S 2323 = S 1313 = S 1 21 2 = C 3 
S 1 1 22 = S 1 133 = S 2233 = C 2 
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where 

C 1 = (Of + C* - 2GC*)/C(C*) 2 + C*q - 2(C|) 2 1 

c 2 = (C*G - C*)/C(C*) 2 + C*C* - 2(C|) 2 ] 

C 3 = {2F 122j 1 C03 + y*/(u' - y*)> 

C* = GF m I 1101 + ( G + 1 ) p i 22 ,1 c0:l + H 

C 1 = F 111 S 1 C0;1 + 2GF 122,1 F0:l " F 
F = -(X* + 2y*)/G 

G = (X* - X*)/C(X' - X*) + 2(y 1 - y*)3 
H = X/G 

Following the theory developed in the previous sections, the 
effective moduli and mass density are found to be 

P* = P + fApD ( 86 ) 


X* - X + fCAX(A nil 2A] 1 22 ) + 2AyA-| i 22 H (87) 

y* = u + fAy(A-j 2 i 2 + A-j 22 i ) (88) 

K* = K + fC(A lin + 2A ]122 )aX 

(89) 

+ (2/3 C3A-J "j 22 A-J212 ^1221^^ 
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It is clearly seen that the velocities are dispersive. At frequency 
range above that of the Rayleigh limit, this phenomenon is pronounced. 
From Figs. 21-23, the bulk moduli, shear moduli and longitudinal veloci- 
ties are shown as functions of volume concentration of spherical 
inclusion materials for the cases of aluminum spheres in germanium, for 
different dimensionless wavenumbers eta. For a given fixed concentra- 
tion, the moduli K*, u* and velocity V[_ and Vj are increased as the 
dimensionless wavenumber ct a is increased. The dispersiveness of effec- 
tive shear modulus is minimal and that of effective bulk modulus is more 
pronounced. Figs. 24,25. 

As an example of application to detect localized damage by void 
nucleation, let all small voids be locally nucleated within a localized 
small region ft of radius R, Fig. 26. The effective moduli of this 
composite can therefore be obtained from Eqs. (71,72). If void nuclea- 
tion outside the region ft can be ingored, then the scattering of the 
composite sphere can easily be obtained. Using the computer program 
developed in [24D, the scattering cross section for a composite sphere 
consisted of small voids in titanium is displayed as a function of 
dimensionless wavenumber for different concentration of voids, Fig. 27. 
It is noted that as the volume fraction of voids insider ft is changed, 
the effective properties, p*, x* and u* are also changed. Hence the 
attenuation effect is pronounced as the concentration of voids is 
increased. The scattering cross section, which is essentially propor- 
tional to the attenuation C193, increases with increasing concentra- 
tion f r . It appears that these curves can be used to locate and 
calibrate porosity in a structural component. Dynamic effective 
properties for this material system are presented in [15]. 
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CLOSING REMARKS 


The mechanics aspects of the characterization of microfracture and 
microdamage by ultrasonics are studied by first looking into the scatter 
of elastic waves by a flat ellipsoidal crack and then by seeking an 
average measure of damage. The work was a part of a three-year program. 

The solution to the direct scattering of a flat ellipsoidal crack 
is presented by using the extended version of Eshelby's method of equiv- 
alent inclusion and a limiting concept. The solution is thus obtained 
by collapsing an ellipsoidal void to a flat crack, say taking a^O. 

The orientation of the crack is assumed to be known. The solution form 
is analytic in incident wave frequency and is in terms of the eigen- 
strains and eigenforces which are governed by the incident wave charac- 
teristics and the equivalence conditions, Eqs. (6,7). The solution 
agrees with the Rayleigh limit [8] and goes beyond it. The solution 
appears to possess a range of validity along the axis of dimensionless 
wavenumber aa^ less than 2 *. 

There are identifiable critical frequencies at which the scattered 
displacement amplitudes become infinite in value. For any given aspect 
ratio of the crack axes, the difference in critical frequencies at sub- 
sequent peaks is inversely proportional to the crack size. A procedure 
for ultrasonic crack sizing is thus suggested and described as follows. 

First, it is assumed that the orientation of the crack plane is 
known or can be determined by finding the direction of maximum scattered 
energy. The differences in frequencies at peak values in the frequency 
spectra at different look angles can then be used to determine the 
aspect ratio and the crack size. The details of the inverse problem of 
crack sizing should be a research program by itself, 

Other areas of research that should be done and can be done are the 
determination of the on-set of microcracking due to orientation and 
geometry by using the crack opening displacement and stress intensity 
factor. Since these quantities can be written in terms of the eigen- 
strains and eigenforces, they can easily be related to the scattered 
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displacements. Earlier references on microcracking in ceramics can be 
found in [16,17,32,37]. Detection and determination of subsurface 
cracks are also of substantial interest in the non-destructive testing 
(NDT) aspect of the science and technology of fracture. The use of the 
concepts in Section III. 2 for possible characterization of transducer 
response is also of interest [33]. 

The velocity and attenuation of ultrasonic waves in two-phase media 
are studied by using a self-consistent averaging scheme. It is required 
that the effective medium to possess the same strain and kinetic energy 
as the physical medium. The concept of volume averaging for physical 
quantities is employed and the solution depend upon the scattering of a 
single inhomogeneity. The thoery is general in nature and can be 
applied to multi -component material system. Since the scattering of an 
ellipsoidal inhomogeneity is known, the average theorem presented in 
this report can be used to study the velocity and attenuation of 
distributed inhomogeneities of shapes such as disks, short fibres, etc. 
The introduction of the orientation of these inhomogeneities besides 
their sizes as in the spherical geometry will necessarily induce 
anisotropy in the effective medium. Fracture toughness and localized 
damage can be studied [4]. 

Results for randomly distributed spherical inclusions of radius "a" 
are presented. Effective moduli and mass density are found to be 
dispersive. The case of a simple model of localized damage is studied. 
Since it is well known that porosity is directly related to the strength 
of rocks and ceramics it appears that the theoretical study of velocity 
and attenuation in two-phase media may be a viable means for data 
analysis in ultrasonic evaluation of dynamic material properties'*’ for 
composite bodies [23,36]. Manufacturing processing, such as rolling, 
sheet metal forming, drawing, etc. often involves plastic flow and 
fracture in the material. Porosity and/or plastic stains induced 


‘'’These are defined as material properties that are obtained by 
using ultrasonics. 
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or contained in the material introduce residual stresses and anisotropy 
in the material and thereby limit the amount of deformation to fracture 
with a directional dependence. Continuous monitor of (1) current global 
moduli, strength and fracture toughness and (2) localized damage such as 
necking may be of importance in the design and optimization of manufac- 
turing procedures. One convenient means for such conti nupus monitor is 
via ultrasonic velocity and attenuation methods. If effective moduli 
and associated phase velocity and attenuation are determined for 
identifiable damage parameters, then the information can be used to 
reconstruct the size and shape of an internal damage zone. Together 
with well developed damage theory, correlation relations with sound 
theoretical basis such as that described in [3,4] will lead to 
prediction of failure or optimum design of processess. 

Other possible application may include soi 1 -structure or fluid- 
structure interaction problems [34,35] where a combination of analysis 
and numerical approach may be involved. The development to cover large 
strain formulation may be needed. 
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Fig. 1 Microcrack process zone 
within the vicinity of 
the crack tip. 




Fig. 2 (a) Transformation + microcracking, 

(b) transformation + microcracking 
+ twinning. (From F. Lange, 1982) 
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Fig. 3 A scanning electron micrograph of a 
typical triple point void in 
sintered A^Og. (From A.G. Evans.) 
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Fig. 4 Inhomogeneous media with multi -components. 
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scattered wave 



Fig. 5 An ellipsoidal inhomogeneity embedded in 
elastic matrix 
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Fig. 6 An elliptical crack: 83 — 0 . 
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Fig. 8 Back scattered data. 
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Fig. 10 Back scattered data. 
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Fig. 12 Back scattered data. 
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Fig. 13 Back scattered data. 
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Fig. 14 Back scattered data. 
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Fig. 16 Back scattered data. 
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Fig. 19 Back scattered data. 
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Fig. 20 Back scattered data. 
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Fig. 25 Aluminum spheres in germanium. 


61 



Fig. 26 A schematic diagram of localized 
damage in an infinite solid. 
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Scattering Cross Section 




APPENDIX A 


Volume Integrals Associated Vlith the 
Inhomogeneous Helmholtz Equation 

Cylindrical Region; Rectangular Region 


I. INTEGRATION OVER A FINITE CYLINDRICAL REGION, Fig. 1 

The integrals, Eqs. (12,13) [1], see also [2], are of either one of the 
following forms: 

4>° = /J7(x') p Cy') q U’) s dv’ 

n 


4> S = ///o(x',y',z’) r r 

n ax' ay' ax' 




(i) 

' sinar ’ , 
1 r’ ; 

dv* 

(2) 

j-cosaY’ i 

dv ? 

(3) 


4> c = ///o(x',y',z') t r — 

n ax' ay ,K ax’ 


Letting x' = x', y'=Ccos0, z'=csin9 and dv'=dx'dy'dz'=Cdi;d0dx' , these 
integrals can be further evaluated as follows: [3] 

(a) 

$° = ///(x') p ( y , ) q ( 2 ') s dv' 
n 


. C°-1)M(»-1)H(2») if all 

(q+s) ! ! ''2+q+s l + p 


even 


(4a) 


if any one of p,q,s is odd. (4b) 


where according to the definition of factorial, 


(q ~ i ) • • = 


(q-l)l 


1 ), 


= 1-3-5. . . (q-1) , q even 


(^) 


! = 2-4-6. . .s(s+2) ... (s+q) , q,s even 


(q-s) ! I = 2 

and a is the radius of the cylinder, C the length. 


(5a) 


(5b.) 
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(b) n=0 , b s : 


» S ■ /|/(x') p (y) < ’(=') s 


m=o 


2m+l - 

= fff(x') p (y') q (z') s l (-l) m ” 1 ) , ( r ’) 2m dv' 

ft 1 ; ‘ 


00 2m*l 

r ( n m a 
u v-7- . i x , S 

m=o ( 2m+ D ! m,p 


where 


S m p = ///C* , ) P (y') q (z , ) S (r') 2n dV 

^ ft 


Using the multinomial formula, [4] 

(r-) 2n = (x ,2 +y ,2 +z ,2 ) m = l — r ; ! |m , (x 1 ) 2m l (y ' ) 2m 2 (s') 2ra 3 

1 2‘ 3‘ 


the integral 


S-- = I 


m! 


ra »P « m m m. inulmj 

m. ,m n 12 3 ft 


///(X ' ) 2m l +p (y • ) 2m 2 + q (;•) 2ra 3 +s dv • 


= 7 gj (u-1) - ! C 1~) ! I (2tt) r a. 2+u+ '\ r 2 ?. 1>X 

m. WV ^v)!! (“ 1 ^) 

1 2 o 

1 1 1 • i • » • , all even 

= 0 if any one of A,y,v is odd 

where (u-1)!!, (v-1)!!, (u + v) ! ! are defined as the same as (3a) and (Sb) , 
and, 

X" = 2m. +d 
1 - 


u = 2m. -<-q 


v = zn.+s 


m = 


( 6 ) 


(7) 


( 8 ) 


( 9 ) 


( 10 ) 
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(c) n=0, $ : 

* c = ///(x') p ( y*) q (x') s S2| Sll d V 

h 

® 2m , n 2m 

= ///(x') p (y ( ) q (x’) s I dv ' 

A m=o 


- Zc-i) 

m=o 


o a 


2m 


(2m) l m,p 


(ID 


where 

C m n = /// Cx') P (y’) q (z , ) S dV' 


ml 


L m,!m 'm_! 
ni l ,m 2 ,m 3 12 3 


/J7(x') X (y') w U’) v pp- 


(12) 


X = 2m^+p 
y = 2m2+q 
v = 2nu+s 


m +m +m„ = m 

1 2 j 


(13) 


The integral on the right hand side of Eq.(12) can be evaluated as 
follows: 


J7J’(x , ) x (y') w (z') v ~r = 

a 

s lu4i!!i^ 

(y+Vj : ! 


= 0 


if u,v all even 


if any one of y , v is odd. 


(14a) 

(14b) 


In the above expression 

tCk-D 


t — 1 ■“ *» 




( 15 ) 
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where 


k = 1+u+v 

is odd , and « , 

c^) rhii) 

fk-l)II = 2 ^ ^ 2 M = 2*4. . .(k-1) 


(k) I ! = . . — = 1*5... (k) 

t^rh±±) 

2 1 2 1 ! 

The integrals 


/ (x , ) X+ 2 t ' 2 (a 2 +x' 2 ) 1 / 2 dx' 


, \ c.,,*!* 1 Ck ...(ki-ZVi) 2t r 2 kl -2 tl .l 3/ j 

k,-l _ 4, 1 J (k 1 +2)k 1 ...(k 1 -2t 1 +4) '■ a - ) (2, - ) (a +i J 


1 * t x =l 


2 (kj"l)!! kj rz T 2 -1 ?. 

+ (”1) (k^+2) (k ) !! 2a ( 4ya *’’ L * a ^ sh J ) if X,y,v are even 


if any one of X,u,v is- odd. (18) 


where 


k L = X + 2t - 2 
is even , and 

Ck.-l) ! 

(k - 1)!; ■ 4 k 


= 1-3-5. 


Ck 1 -1) 


(20a) 


— — ^ l 

(kj]!! * 2 2 = 2-4... (kj) 


is even 


The integral 


/ (xM^ka 2 *^ 2 ) 1 ^' 


= / (X 1 ) 2 (a 2 + x' 2 ) 1 / 2 dx' 
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^2 

“ 


k 2 *fi 


? r ^2^^ ^2" ^ 2t^-2 k--2t-+l ,/^ 

^ ( -» - rv^ 2 ...(k 2 -2t 2 V (a) ' •“> • <* 2 *‘V 2 




2 (k-,-1)!! ^2 n, r ^ „io 

<-» (k 2 .2?-( k pil 2 U /a^ 2 * a~sh |) ifX.u, 


v axe even (21a) 


= 0 


if any one of X,u,v is odd. 


(21b) 


where 


k 2 =X+k-l=X+u+v 


(22) 


The integral 

?* , ,.X+k, . 

J (x 1 ) dx* 

-£ 


22 


X+k+1 


X+k+1 


= 0 


if X,u,v are even (23a) 

if any one of X,u,v is odd. (23b) 
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II. INTEGRATION OVER A RECTANGULAR PARALLELEPIPED, Fig. 2 


In the case of a rectangular parallelepiped, the integrals in (14) 
(15), (16) can be evaluated as follows: 

(a) 

* 0 = ///U') p (y , ) q ( 2 ') s dv- 

a 

- (p.i) (q?D (s.i) - M q * 1 (b) s * 1 m p * 1 if p.q.= •»«* 

= 0 if any one of p,q,s is odd. 


where 2,,a > b is length of the rectangular parallelepiped toward x' f y',z f - 


direction, respectively, 
(b) n=0, $ s 


* s = J7/(x') p (y') q (z') s dv' 

n 


2m+l - 

m a 2m 


= J/J(x*) P (y') q (z') s l (-i) m “— r (r) 

m-r ' 


nr=o 


2m+l 


" X C_1) (2m+l) ! S m,j 


m=o 


and 


S 

m,p 


///(x') p (y') q (2') S ( r') 2m dv' 

n 


. 2m, +p 2m,+£ 2m +s 

l izitirrW*'’ ‘ (I,) dv ' 


“I 4 2 3‘ 0 

1 2 o 

„ . , , ji+1 , v+1 X+l 

r 8m! (a) (b) (&) 

_ “ _ m.Jm'm _! (X+l) (y+1) (v+1) 

yl&2 


Hi— 0 y 1 y 2 | a a • 
= 0 


if X,u,v all are even 
if any one of X,u,v is odd. 


(24a) 

(24b) 


(25) 


(26a) 

(26b) 


where 
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m = m +m +nu 

1 2 j 


X = 2m^+p 
U = 2m 7 +q 


v = 2m_+s 
o 


( 27 ) 


(c) n=0, <t> . 

* ///(x') p Cy') <1 (z , ) s -£2^£l dv’ 

r\ r 


2m 


= ///(x') p (y-) q (z') s 


2m 


= l (-D 


m a 


2m 


(2m)! m,p 


(28) 


and 


m. 


P’ m I „ r i S iU i/W(r-)*U')”? 


£ (x ' )X I/ <>">“</ dz')dy]dx' 


m. ,m 7 ,m- m l- w 2- ,u 3* -4 -a -b /2“^ 
i 2 o /ct +z 


(29) 


where 


X = 2m^+p 
u = 2m 2 +q 
v = 2m 3 +s 
a 2 = x ' 2 +y ' 2 


= m 

1 2 j 


(30) 


The integrals in (34) can be evaluated as follows: 


(i) b 

/ 


(S') 


-b /~2 “ 

/a +z ? 


dz ’ 


_ _ r , v v-2t+l c .2 f 2x L2 72 72 f ,.v/2 2(v-l) ! * f 2 ,2.v/2 ,-l b 
= 11 v ( b ) (X' + y' )/b +x' +y ' + (-1) — ( x + y 1 sh T=t 


/x-^v^ 
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=0 


if v is even (31a) 

if v is odd. (31b) 


where is defined as an operator as follows: 

n _ _2_ y /2 f n t+l (v+l)(v-l)...(v-2t+3) 
v ■ v+1 ‘j *•" J v(v-2)(v-4)...(v-2t+2) 


(32a) 


(v- 1) ! ! = 


(v-l)l 


= 1*3. . . (v-1) 


(v)!! 


! - 2^ ! . 2- 


2*4. . .v , v is even 


(32b) 


(32c) 


(ii) 


/* /V)"u') u . 

” a -b /x , “+y ,2 +z' 2 


« n v I c-d v/2 r7rr (b)V ‘ 2t * lcx ’ )2ni V*’) 


n l’ n 2 


12 


* : I, c-D v/2 ipijr 2( cv)i!^ b)(x ' )211 V*' 5 


1 l ,1 2 


v/2_ij_ [.•l)!! w i 1 , l) ii l !.<*■) 


y f_n v > ^ i v — lll 

. L. 1 iJ 3(v) ! 1 

H’ 1 ! 1 2 

CO 

I I (-1) 


v/2+j i! 2(v-l) ! ! (2j) ! 


i l ,i 2 j =2 


W (v)!! 2 2j (j!) 2 (2j+l) 


(b) 2;i+1 (x') 211 I 4 (x') (33) 


where 


3. 2n2 + h /*5 5 2* 

[ 1 (x') = / (y») Ar+x'^y'* dy' 

-a 




• (a>/a 2 +b 2 +x' 2 + (b 2 *x ,2 )sh _1 a =) 

/b 2 +x’ 2 
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= 0 


if u is even 
if u is odd. 


In the above expression 

n - 2 Z ( 2 r n Z l + 1 (5 + l)(5-l)...C5 + 2t 1 +3) 

C " v+1 t ^i (5+2)5. ..(C-2t 1 + 4) 


(34a) 

(34b) 


(35) 


(5-1)!! = 1*3. ..(5-1) 5 is even 
(5+2)!! = 2-4... ( 5 + 2 ) 5 is even 
5 = 2n 2 +u 


(36 a) 
(36b) 
(37) 


The integral 


a 2 i-+u 

i,(x') = / ( y’) 


dv f 


/x f 2 +V ' 


= n, 


- *3 1 ]~2 J ^]/ 2 2(C^-1)!! 2 ^]/ 2 _ ^ a 

. (a) 1 J (x' 2 ) + (-1) -±— (x^ 2 ) sh 1 f- 

’ i U 1 J, ‘ 


= 0 


if u is even (38a) 

if y is odd. (58b) 


where U is defined as the same as (52a) , only if replace v with £ ^ and 


5i = 2 i 2 +u 


(39 a) 

(5 X -1) ! ! = 1-3. ..(5 r l) 

is even 

(39b) 

(spn = 2-4. ..fcp 

is even 

(39c) 


The integral 


I 3 (X') = / 


(/’) 


-a 


*{x f 2 +y 1 2 ) 0 

(X- 2 ) 


= dy' 


^l"“ t 4 + ^ 2 t 4"^//2 2 

= n^(a) 1 (x' Z ) //a +x' Z 

^I 2 


1 1 

T~“2 (5 ,-1) ! ! 2 2 1 -la a 

+ (5tW (X } Csh tr-^=7=) 


I^^ + X' 
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x 


if u is even 
if y is odd. 


wnere 


*1 


2 5 1 ; 1 «,*• 

v 1 vi £ - l) 


CCj+DCSj-l) ...CC r 2t 4+ 3) 
(? 1 -2)(? 1 -4)...(C 1 -2t 4 ) 


(40a) 

(40b) 


(41) 


The integral 


I 4 (x') = / 


(ill 


5l 


-a /(x '-+ y ' 4) 4j 


dy' 


Ci“2u- + 1 « ^7” 1 


= n" (a) 1 2 (x’ 2 ) 2 //(a 2 +x' 2 ) 2;i ' 1 

^1 

i” 1 ^ ^V 2 2(C r l)!!(2) to 

+ t ( ' 1} (S 1 -2j)« 1 -2j-2)...(2-2j) (2j-l)(2j-'3T::.(2j-2t *1) 

° 5 1 . 

• (a)(x- 2 ) 2 • t °' 1 Lai^ 2 ) C ° 


(j-1) (j-2) . . . (i- t o) 


(- 1 ) 


/ (a^+x* 2 ) 2 3” i 

V 2 2C€ r l)N 


C5 1 -2j)(5 1 -2j-2)...(2-2j)C2j-l) 


2 2 

■(a) (x'V 


-- 1 


/ (a 2 +x r2 ) 2 ^ 


= 0 


if y is even 
if y is odd 


where 

7 ? l/2 t,+l (c,*l) (c,-l) • • • (5,-2t-+3) 

n" = — v r_n * t 1 £ 

? 1 V 1 t -1 (C 1 -2j)(5 1 -2j-2)...(C 1 -2j-2t 2 *2) 


C 42a) 
(42b) 


(43) 
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(iii) 


//J(X-)V) U ( Z') V — dV> : 

ft /x ,2 +y ,2 +z' 2 


= £ (-l) v/2 ? - ! , (a) g 2tl •(b) v ‘ 2t+1 -I. 

v5 ni ,n 2 n i !n 2 ! ] 




J. (-» V ' 


v/2 i! 2Q-1) 


! ! , , ? l‘ 2t 2 +1 


1 i. ,i 


V ( _i) 2 _jl 2Cv-_ l)l! 2(? r 1)!! (b) tI 
, S L J W 0)N (5J!! l J X 4 


. n . y f n v/2 i! (v-Dn ^A^V 1 ^^ 

V . * ( 1} i.li,! 3(v)ll Ca) Cb) X S 

1 » 1 2 1 ^ 


5,+v-l 

p “^2 i- 0-1)" 2 ^ 1 -1)!! 3 

W 3 w |! ‘v 2 > !l tb) '“ 6 


K.l. U-» rfirr 


2 i! 2(v- 1) ! ! C 2 j ) ! 


5,-2t,+l 


'! i 1 »i 2 2 =2 


(v)!! 2 2j Cj!) 2 (2j + l) 


4 j4 


y y 'j ( 4^2" i! Ihzllll l-i’— 

, l , (v)!! (5,-2j) CC,-2j-2) . . . (2-2j) 

1’ 2 - 1 o 


(2j)! (j-l)(j-2)...(j-t 0 ) 

2 2j (j!) 2 (2j + l) C2j-l)C2j-5)...C2j-2 V l) 


■(a)(br J -I. 


. v 5l 

J 2' 2 i! 20-1) ! 1 (2 j ) ! 


y y r n 2 2 l! IQ) ; 

i^J^ ’ W 2 2 ^ (j ! ) 2 (2j + l) 

. 2(? 1' 1)!! (aHb4 +1 -I 

(5 1 -2j)(C 1 -2j-2) . . . (2-2j) ( 2 j - 1 ) w lDJ 9 


The integrals in (44) can be evaluated as in Appendix B. 
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appendix b 


The I-Integrals and it-Functions 


(a) The I -integrals 

(i) 

h = J 

•Z 


1 1 = j Cx’) X+2:i Cb 2 +x' 2 ) 1 1 /(a 2 *b 2 +x' 2 ) 3 dx’ 


2 r. 


Vi t ^ 1 n-2t T +l — r~ 

- * I r ,-r~, (b) 1 -(a 2 *b“) -U) /Ca 2 + b 2 * 2 ) S 

1 1’ 1 2 1 2 


* I r (*1) n/ *" fTT^T ■^73j7T^*(b) *"* 1 • (a 2 +b 2 ) n//2 [iz- /(a 2 +b 2 +Z 2 ) 3 
r i’ r 2 1 2 


t 7 7 . *■» 9 7 7 i o 

+ -rCa^+b 4 ") (z/a^*b*W~ + (a**+b^)sn" 


/a 2 -b 2 


0] 


ir \ is even 

=0 if \ is odd 

in which the t- function is defined in fbj of the Appendix. 

n 

n-n i = i* + r» 


X+2n, , 


I! = ' a(x') 1 (b 2 *x' 2 )’ /2 -/a : *b 2 .x'" dx’ 


-i 


j. t 26 | t 1 -l / -7 9 - 

= H 7 a » (a) (b) * • (a~-b~ j " *(£) ** * * (a“*b 

n i a “a ~ 1 * ’ 7 * 


“1 A a 

x - ^ I w 9 


n j/2 , 2(n,-l)!: 

Y /* -I A x J * 


C-i) 


8 1 >5 2 


e i :3 2 : < n r 2)!I 


2S 1 7 7 r ‘]/ 2 

TTTr-CaJCb) (a~+0 


(z/a 2 +b~+ 


7 ? -i i 

+ (aMT)sh 


= 0 


r'a^b 4 - 
if X is even 
if X is odd. 
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I X+2n, 

1 2 = / Cx f ) 1 (b 2 +x’ 2 ) 35/2 sh' 1 a Hr' 

-1 / b 2 +x <2 


n -/2 

p , y t. + i 
TrrC i (-D 


20, n.+l ru/^ 

Z £• a(b) V) J _ J. p e! , J 7 

n.+l 5 ^ p i e i ^ “ 

1212 o e'e,! e l* 2‘ t =1 

1 ^ 6 

f 7 n_-2t,+i 

jo> 2 ) 6 (&) 


n--2t-+l 


+ (-1) J (b) 3 -t - 1 A )(a) 3 (b) -'l 


V p j l C2 ^P ! 

l l (-D — 


e! ,.2il +l 


S l’ e 2 •’l -2 


22jl O' 1 JI(2j 1 +lJ e i !e 2 !(a) 


2 1 r l-“ 2 ty+l 

— <-» 3 

(b 2 +£ 2 ) 1 


^3-1)!! 


(n 3 -2j 1 *l)(n--2j-l)...(3-2j 1 ) 


n 3 JV 1 ? i o 2e, 

* (b) 3 C!T (i) + (-I) | tg -1 |) ] (b) 1 


if X is even 


if X is odd. 


(iii) * A+2i 7 +2t,-2 — — 

I? = / (x') 1 3 /a 2 +x' 2 dx' 


2 t 5 ' 1 r| 2 " 2t s +1 , 3 — n,/2 2(n,-l)!! _ r> n 

) CO /(a +n ♦ (- 1 ) _L__ {a 2 ) 


= n (a ) 00 

n 2 


(Jl/a 2 + Jl 2 + a 2 sh 1 -) 


if X is even 


if X is odd. 


I X+2i_+C 1 

I 4 = / (X') W 1 fr dx' 

-i 

a 2 t 15 -1 V 2t lS + 1 n g /2 2(n g -l)!! n g . 

■ ^ty* > w /a_ * 1 * ( ' I) -yrr w sh 


if X is even 


if X is odd. 
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oj | Oj 


(V) 


l A+Zi.+Zt.-Z 

; ex’) 1 4 

-l 


dx' 

/a^+x'^ 


• n (a) 
n 8 


h^ 1 


U) 


n 8- 2t 14“ 1 


2.,Z 


• / a" + i 


n s /2 2( V 1,U „ s /2 , 

( - 1} ~ (n 8 ) H (a > sh X 


= 0 


if X is even 
if X is odd. 


Cvi) 


I = I» - T" 
l 6 l 6 l 6 


% X+2i 1 +C 1 -2 . 3. i | 

I' = / ex') 1 1 sh' 1 3F ^ 

b -l 


» pi u 2 )' 16 ' 1 -^)" 16 ' 2 ' 16 * 1 -^ 2 * (-D 


n 10* 11 "10 


n 10 /2 2 ( n 10‘ 1)!l 

(n 10 )l! 


(a) “ah' 2 i] 


if X is even 
if A is odd. 


I” = / (x>) 


l X + 2i^+£^-2 


-£ 


/a^+x' ^ 


dx' 


n r -2t Q -i 


H (a) 9 -C*) 5 9 -/a 2 *l 2 + (-1) 

"e 


V 2 2Cv 1 ) ,, ,.,V 2 . u -i 

(n s ) ! ! C 5 h 


= 0 


if X is even 
if X is odd. 


. . . . £ X+2i. +2t~-2 , t 

Cvil) I_ » / (x-) 1 2 dx -- 

-l /(a2+x'2) 2 j _1 


_ t in -l m V 2t 10 +1 n fi /Z 

= n (a 2 ) 10 ■ <& . — - + (.i) 6 


Z(n 6 -1)!! 


/(a^^) z j-3 


Cn 6 -Zj-2)(n 6 -2j-4)...(-2j) 
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j ' 2 3 kll (j-2)(j-3)...(j-k n -l) (a 2 + 2 ) kl1 

C2i-3)a 2 /Ca 2 +i 2 )^‘ ,i ^ k,%l ^ ' 5) (2j ' 7) * ‘ ’ ^ ' 2k H' j) ,,_2, k ll ] 


a 6 i 


nr 


(4a fc ) 


= 0 


if X is even 
if X is odd. 


(viii) Z A + 2i 1+ 5 r 2t 0 -2 (a 2 +x ,2 ) o^, 

I = / (x') 


-% 


* / (a 2 +x ,2 ) 2 J‘ i 


= n 


. , 2n t 12" 1+ ^ l,., n 7~ 2t 12** 

4! (a") U) 


4! 


n 7 . $,!*,! 

7 9 1 >'t i 2 1 ~ 


/ (a 2 +H 2 ) 2 J " 2 


L 4 . 1 o _ i 
*,,*■> 2 


(-D 


n ? /2 2(n 7 -l)!! 


(n 7 -2j-2)(n 7 -2j-4)...(-2jJ 


(a 2 )' 7 


n 7 /2+4 . 


j : 2 8 kl3 (j-2) ( j - 3) . . . (j-k 13 -l) (a 2 + i 2 ) kl3 


* a + — 

(2j -3) a 2 /(a 2 +2. 2 ) 2 J + k.t-l (2j-5) (2j-7) . . . (2j-2k -3) (4a 2 ) k 13 


if X is even 
if X is odd. 


dx' 


(ix) l X+2i. +£. -2 

r g = / (X') ; = = -i, = 

-l /(a 2 +x' 2 ) 


2 V 1 V 2t 12 + 1 

n Ca 2 ) X - U) 
n 8 /(a 2 +£ 2 ) 2 J -1 


n 8 /2 2(n 8 -l)J! 

(n 8 -2j-2)(n 8 -2j-4)...t-23j 


a \ 


(2j-3)a 2 /(a 2 +Z 2 ) ^ "~' 5 


(1 + 


j " 2 8 kl2 (j-2)(j-3)...Cj-k 12 ‘ 1 )(a 2 +r) 12 

k 12 =l (2j -5) (2j-7) . . . (2j-2k 12 -3) (4a 2 ) k l 2 


= 0 


if X is even 
if X is odd. 
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(b) The tt- F unctions 


The functions in the I-integrals listed in (a) of this Appendix are 
defined as follows: 

(i) 2 n ^ 2 (n+1) (n-1) . . . (n-2t-+5) 

n n ” ri+r ^ (n+4) Cn+2)...Cn-2t-+6) 

5 J 

n = A. + 2n^ + 2T2 (is even) 

T + r 2 * F = V 1 


(ii) n as the sane as (35) , if reDlace 5 with n, 
'i * i 


= \ + 2n^ +29 ? (is even) 


Sj + e 2 = 9 = 5/2 


(iii) 


n-/2 


n_+l 
3 t-=l 


l (-D 


t 7 +1 (n-+l) (n -- 1) . . . (n--2t_+3) 


(n 3 -2j 1+ i) (n.-2j r l) . . . (n 3 -2j r 2c ? +3) 


4 Jl " r \ „V1 _yl [3-2U 1 +l)][3-2j 1 j...[3-2(j 1 -t s +2)] 
n 3 5 - 2 ^l +1) t a =l U 1 -D(j 1 -2)...(j 1 -t 8 ) 


n- = A + 2n^ +• 2e 9 (is even) 


e l + e 2 e " 2 5 


(iv) II as the same as (35), if replace 5 with 

n 2 2 


n 2 = ^ + + 2t, “ 2 (i s even) 


(v) IT as the same as (32a) , if replace v with n Q 
n 9 * y 


rig = A + 2i^ + (is even) 
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(vi) 2^ as the same as (52a), if replace v with n q 

3 * 6 

n s = X + 2i^ +■ 2t 4 -2 (is even) 

(vii) 2^ as the same as (52a) , if replace v with 

= X + 2i^ + ^ - 2 (is even) 

(viii) H as the same as (32a), if replace v with 

n 5 0 

rig = X + 2i^ * ^ -2 (is even) 

(ix) , rl ^{" t 10 *l (n 6 *l)(n e -l)...(n 6 -2t 10 +3) 

\ ' V r tj.l C ' 1) Cn 6 -2j-2)Cn 6 -2J.4)...(„ 6 -2j-2t 10 ) 

n 6 - A + 2i^ +2t 0 - 2 (is even) 

(x) n as the same as II , if replace n* with ru 

n_ n ’ o / 

/ 6 

ru = X + 2i ^ ^ +2p ? - 2t Q - 2 (is even) 

+ *2 = * = t Q 

(xi) a i as the same as , if replace with rig 

8 6 

ng = X + 2i, + Cj - 2 (is even) 


83 









National Aeronautics and third-class bulk rate 

Space Administration 

Washington, D.C. 

20546 


Official Business 

Penalty for Private Use, $300 


Postage and Fees Paid 
National Aeronautics and 
Space Administration 
NASA-451 





POSTMASTER: 


If Undeliverable (Section 158 
Postal Manual) Do Not Return 



